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Abstract — We present a mathematical model for communica- 
tion subject to botli network interference and noise. We introduce 
a framework where the interferers are scattered according to 
a spatial Poisson process, and are operating asynchronously in 
a wireless environment subject to path loss, shadowing, and 
multipath fading. We consider both cases of slow and fast- 
varying interferer positions. The paper is comprised of two 
separate parts. In Part I, we determine the distribution of 
the aggregate network interference at the output of a linear 
receiver. We characterize the error performance of the link, in 
terms of average and outage probabilities. The proposed model 
is valid for any linear modulation scheme (e.g., Af-ary phase 
shift keying or A/-ary quadrature amplitude modulation), and 
captures all the essential physical parameters that affect network 
interference. Our work generalizes the conventional analysis of 
communication in the presence of additive white Gaussian noise 
and fast fading, allowing the traditional results to be extended 
to include the effect of network interference. In Part II of the 
paper, we derive the capacity of the link when subject to network 
interference and noise, and characterize the spectrum of the 
aggregate interference. 

Index Terms — Spatial distribution, Poisson field, aggregate 
network interference, error probability, stable laws. 

I. Introduction 

IN a wireless network composed of many spatially scattered 
nodes, there are several fundamental impairments that con- 
strain the communication between nodes, including thermal 
noise and network interference. Thermal noise is introduced 
by the receiver electronics and is usually modeled as additive 
white Gaussian noise (AWGN), which constitutes a good ap- 
proximation in most cases. Interference, on the other hand, is 
due to signals radiated by other transmitters, which undesirably 
affect receiver nodes in the same or in a different network. For 
simplicity, interference is typically approximated by AWGN 
with some given power f\\. However, this elementary model 
does not completely capture the physical parameters that 
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affect interference, namely: 1) the spatial distribution of nodes 
scattered in the network; 2) the transmission characteristics of 
nodes, such as modulation, power, and synchronization; and 
3) the propagation characteristics of the medium, such as path 
loss, shadowing, and multipath fading. If, instead, a spatial 
Poisson process is used to model the user positions, then all 
these parameters are easily accounted for, and appear explicitly 
in the resulting performance expressions. 

The application of the Poisson field model to cellular 
networks was investigated in Q and later advanced in Q. 
However, these papers either ignore random propagation ef- 
fects (such as shadowing and multipath fading), or restrict the 
analysis of error probability in non-coherent FSK modulations. 
In other related work ||4l, it is assumed that the different 
interferers are synchronized at the symbol or slot level, which 
may be unrealistic in most situations. In ||5l, the authors choose 
a different approach and restrict the node locations to a disk or 
ring in the two-dimensional plane. Although this ensures that 
the number of interferers is finite, it complicates the analysis 
and does not provide useful insights into the effects of network 
interference. In l|6J-|i8J, the authors analyze coexistence issues 
in wireless networks, but consider only a small, fixed number 
of interferers. Lastly, none of the mentioned studies attempts 
a system characterization that incorporates various metrics 
such as error probability, channel capacity, and power spectral 
density. 

In this two-part paper, we introduce a more realistic frame- 
work where the interferers are scattered according to a spa- 
tial Poisson process, and are operating asynchronously in a 
wireless environment subject to path loss, shadowing, and 
multipath fading ll9l- lfT2l . We specifically address the cases 
of slow and fast- varying interferer positions. In Part I of the 
paper, we determine the statistical distribution of the aggre- 
gate network interference at the output of a linear receiver, 
located anywhere in the two-dimensional plane. We provide 
expressions for the error performance of the link (in terms 
of average and outage probabilities), which are valid for any 
linear modulation scheme. We then quantify these metrics as 
a function of various important system parameters, such as the 
signal-to-noise ratio (SNR), interference-to-noise ratio (INR), 
path loss exponent, and spatial density of the interferers. Our 
analysis clearly shows how the system performance depends 
on these parameters, thereby providing insights that may be 
of value to the network designer In Part II of the paper ifTSll . 
we derive the capacity of the link when subject to network 
interference and noise, and characterize the spectrum of the 
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aggregate interference. 

This paper is organized as follows. Section HIl describes the 
system model. Section |III] derives the representation and dis- 
tribution of the aggregate interference. Section |IV] analyzes the 
error performance of the system, and gives plots to illustrate 
its dependence on important network parameters. Section [V] 
concludes the paper and summarizes important findings. 

II. System Model 

A. Spatial Distribution of the Nodes 

We model the spatial distribution of the nodes according to 
a homogeneous Poisson point process in the two-dimensional 
infinite plane. Typically, the terminal positions are unknown 
to the network designer a priori, so we may as well treat them 
as completely random according to a spatial Poisson processQ 
Then, the probability P{n in TZ] of n nodes being inside a 
region TZ (not necessarily connected) depends only on the total 
area A of the region, and is given by lfT4l 



"{n in 7^} 



: — e 



n>0, 



where A is the (constant) spatial density of interfering nodes, 
in nodes per unit area. We define the interfering nodes to 
be the set of terminals which are transmitting within the 
frequency band of interest, during the time interval of interest, 
and hence are effectively contributing to the total interference. 
Then, irrespective of the network topology (e.g., point-to- 
point or broadcast) or multiple-access technique (e.g., time 
or frequency hopping), the above model depends only on the 
density A of interfering nodes 

The proposed spatial model is depicted in Fig. [T] For 
analytical purposes, we assume there is a probe link composed 
of two nodes; the probe receiver, located at the origin of the 
two-dimensional plane (without loss of generality), and the 
probe transmitter (node i = 0), deterministically located at a 
distance tq from the origin. All the other nodes {i = \ . . . oo) 
are interfering nodes, whose random distances to the origin 
are denoted by {Ri\°^i, where Ri < R2 < ■ ■ ■■ Our goal is 
then to determine the effect of the interfering nodes on the 
probe link. 

B. Transmission Characteristics of the Nodes 

To account for the transmission characteristics of users, we 
consider that all interfering nodes employ the same linear mod- 
ulation scheme, such as A/-ary phase shift keying (AZ-PSK) 
or A/-ary quadrature amplitude modulation (7\/-QAM), with 
symbol period T. Furthermore, they all transmit at the same 
power P - a plausible constraint when power control is too 
complex to implement (e.g., decentralized ad-hoc networks). 
For generality, however, we allow the probe transmitter to 
employ an arbitrary linear modulation and arbitrary power Pq, 
not necessarily equal to those used by the interfering nodes. 

'The spatial Poisson process is a natural choice in such situation because, 
given that a node is inside a region TZ, the probability density function (p.d.f.) 
of its position is conditionally uniform over 7?.. 

^Time and frequency hopping can be easily accommodated in this model, 
using the splitting property of Poisson processes 1151 to obtain the effective 
density of nodes that contribute to the interference. 



O Probe transmitter node 
o Probe receiver node 
• Interfering node 




Figure 1 . Poisson field model for the spatial distribution of nodes. Without 
loss of generality, we assume the origin of the coordinate system coincides 
with the probe receiver. 




symbols: 
node i 




D, 



Figure 2. Asynchronism between different transmitting nodes. In the 
observation interval [0, T], a change in constellation symbol of node i occurs 
at random time t = Di, from aiC^^' to a'^e^^i, where a and 9 denote the 
transmitted symbol amplitude and phase, respectively. The distribution of Di 
is assumed to be U{0,T). Therefore, node initiates symbol transmissions 
at times nT by convention, while node i initiates symbol transmissions at 
times nT + Di. 



We do not assume synchronization among interfering nodes, 
but instead consider asynchronous transmissions where dif- 
ferent terminals are allowed to operate independently. As 
depicted in Fig. |2] node i transmits with a random delay Di 
relative to node 0, where Di ^ U{0,T)^ The probe receiver 
employs a conventional linear detectorO Typically, parameters 
such as the spatial density of interferers and the propagation 
characteristics of the medium (e.g., shadowing and path loss 
parameters) are unknown to the receiver This lack of infor- 
mation about the interference, together with constraints on 
receiver complexity, justify the use of a simple linear detector, 
which is optimal when only AWGN is present. 

C. Propagation Characteristics of the Medium 

To account for the propagation characteristics of the envi- 
ronment, we consider that the median of the signal amplitude 

'We use U{a, b) to denote a real uniform distribution in the interval [a, b]. 
''Note that the other receiver nodes are not relevant for the analysis, since 
they do not cause interference. 
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decays with the distance r according to k/r^, for some given 
constant k. The ampHtude loss exponent b is environment- 
dependent, and can approximately range from 0.8 (e.g., hall- 
ways inside buildings) to 4 (e.g., dense urban environments), 
where 6 = 1 corresponds to free space propagation ml The 
use of such a decay law also ensures that interferers located 
far away have a negligible contribution to the total interference 
observed at the probe receiver, thus making the infinite-plane 
model reasonable. 

To capture the shadowing effect, we model the channel 
amplitude gain as a log-normal random variable (r.v.) 
such that S = fie'^^, where G ^ 7V(0, 1)^ fi = k/r'' is the 
median of S, and a is the shadowing coefficient^ Thus, 
the shadowing is responsible for random fluctuations of the 
channel gain around the median path gain k/r^. The multipath 
effect is modeled as fast fading, which is superimposed on the 
path loss and shadowing. Specifically, the fading affects the 
received signal by introducing a random phase (j) ^U{0, 27r), 
as well as an amplitude factor a with arbitrary distribution 
and normalized to have unit power gain, i.e.. 
Because of its fast nature, the fading is always averaged out in 
this paper, both when determining the interference distribution 
and the error probability. 

In what follows, we consider the shadowing (and similarly 
for the fading) to be independent for different nodes i, and 
approximately constant during at least one symbol interval. 
Additionally, the probe receiver can perfectly estimate the 
shadowing and fading affecting its own link, hence ensuring 
that coherent demodulation of the desired signal is possible. 

D. Mobility and Session Lifetime of the Interferers 

Typically, the time variation of the distances {Ri\°^i of 
the interferers is highly coupled with that of the shadow- 
ing {Gi]°^i affecting those nodes. This is because the shad- 
owing is itself associated with the movement of the nodes near 
large blocking objects. Thus, we introduce the notation V to 
denote "a particular realization of the distances {Ri\°°^i and 
shadowing {Gi}°^i of the interferers," or more succinctly, 
"the position of the interferers" In this paper, we analyze the 
following two scenarios, which differ in the speed of variation 
of P: 

I) Slow-varying V: During the interval of interest (e.g., a 
symbol or packet time), the distance Ri of each in- 
terferer is approximately constant, Ri{t) w Ri. Fur- 
thermore, the interferers have a long session lifetime, 

^Note that the amplitude loss exponent is h, while the corresponding power 
loss exponent is 2b. 

*We use Af{p., cr^) to denote a real Gaussian distribution with mean p and 
variance cr^. 

'This model for combined path loss and log-normal shadowing can be 
expressed in logarithmic form 1161 , 1171 . such that the channel loss in dB 
is given by Ljb = ko + ki log^Q r + cr^^G, where G ~ A^(0, 1). The 
environment-dependent parameters (fcoifciifdB) can be related to {k,b,a) 
as follows: fco = — 201og]^Qfe, k\ = 20b, and o-jb = "nTTo"'- The 
parameter (tjb is the standard deviation of the channel loss in dB (or, 
equivalently, of the received SNR in dB), and typically ranges from 6 to 
12. 

*We use E{-} and V{-} to denote the expectation and variance operators, 
respectively. 



transmitting continuously over many symbols. In this 
quasi-static scenario, V varies slowly with time, and thus 
it is insightful to condition the interference analysis on 
a given realization of V. As we shall see, this naturally 
leads to the derivation of the error outage probability of 
the probe link, which in this case is a more meaningful 
metric than the error probability averaged over V [[TSl - 

2) Fast-varying V: As in the previous case, Ri{t) « Ri 
during the interval of interest. However, the interferers 
have a short session lifetime, where each node periodi- 
cally becomes active, transmits a burst of symbols, and 
then turns off (e.g., in a sensor or a packet network). 
Then, the set of interfering nodes (the set of nodes 
that are transmitting and contributing to the interference) 
changes often, and so does their effective position V, 
which experiences a variation analogous to that of a 
block fading model. In this dynamic scenario, it is 
insightful to average the interference analysis over all 
possible realizations of V, which naturally leads to the 
derivation of the average error probability of the probe 
Unk. 

III. Interference Representation and Distribution 
A. Complex Baseband Representation of the Interference 

Under the system model described in Section HI] the ag- 
gregate signal Z{t) at the probe receiver can be written for 

< t < r as 



Z{t) 



kaoe' 



<yGo 



-ao cos(27r/et + ^o) + Yit) + Wit), 



where the first right-hand term is the desired signal from the 
transmitter probe node, Y{t) is the aggregate interference with 



kc 



R 



—Qi cos(27r/ci + + (j)t)u{Di - t) 



+ \j -a[ cos(27r/et + e[ + (j),)u{t - A) 



< i < T, 



and W{t) is the AWGN with two-sided power spectral den- 
sity No/2, and independent of Y{t). In the above equations, 
we use the following the notation: T is the symbol period; /c 
is the carrier frequency; Uie^^' and a^e^^* are rv.'s denoting 
successive constellation symbols transmitted by the node i 
during the interval of interest [0,r] (see Fig. |2]i; and u{t) is 
the unit step function. The overall effect of the path loss, log- 
normal shadowing, and fading on node i is captured by the 
ampHtude factor kaie"'^' / R\, where --7V(0, 1), and by 
the uniform phase We consider that r.v.'s a,;, 0i, G;, Ri, 
aie^^\ a'iC^^i, and Di are mutually independent for a given 
node i, and that the sequences {ai}, {Gi], {oie-'^'}, 

{a^e^^'}, and {Di} are independent identically distributed 
(i.i.d.) in i. 

'since we assume the probe receiver perfectly estimates the phase (po of 
the multipath fading affecting its own link, we can set 0o = without loss 
of generality. 
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The probe receiver demodulates the desired signal from 
the aggregate signal Z{t), using a conventional lin- 
ear detector. This can be achieved by projecting Z{t) 
onto the orthonormal set = ■\/2/T cos(27r/ct), 

-02 (i) = — ■\/2/r siii(27r/cf ) | . Defining the in-phase and 

quadrature (IQ) components Z„ = Z{t)il;nit)dt, n = 1,2, 
we can write 



Zi = 



-aoCos6lo + Yi +Wi 



kaQC 



-ao sin 6*0 + I2 + W2, 



(1) 
(2) 



where Wi and are A/^(0, Nq/2) and mutually independent. 
After some algebra (Appendix IaIi. Yi and I2 can be expressed 
as 



1,2, (3) 



where 



fca, [ai^ cos{e.i + (f)i) + a- (l - -^f) cos(6'- + 0,;)] 



i,2- 



,^ sin(0, + 0,) + a^l - sin{9[ + 0,)] 



(4) 

I- 
(5) 

Using complex baseband notation^ equations (|T|)-(|5l) can be 
further simplified as 



-floe- 



■J a a 



w 



E 



where 



D, 



-a,e>"^ + 1 



and the distribution of W is given b\Fl 
W^A/-e(0,iVo). 



(6) 
(7) 

(8) 
(9) 



Since different interferers i transmit asynchronously and inde- 
pendently, the rv.'s {X^}^]^ are i.i.d. 

The distribution of the aggregate interference Y plays 
an important role in the evaluation of the error probabiUty 
of the probe link. In what follows, we characterize such 
distribution in two important scenarios: the P -conditioned and 
unconditional cases. 



B. V -conditioned Interference Distribution 

To derive the P-conditioned distribution of the aggregate 
interference Y in ©-([SI), we start with the results given in 
). This work shows that in the case of Rayleigh fading, an 



'"Boldface letters are used to denote complex quantities; for example, 
Z = Zi+jZ2. 

"We use A/'c(0, o"^) to denote a circularly symmetiic (CS) complex Gaus- 
sian distribution, where the real and imaginary parts are i.i.d. Af{0, a'^ /2). 



expression of the form of (|8]l can be well approximated by a 
circularly symmetric (CS) complex Gaussian rv., such that 



■NA0,2Vx), Vx=Y{X,,n}. 



(10) 



In II22I . the validity of this approximation is justified both 
by analyzing the Kullback-Leibler divergence and comparing 
the error probabilities in the exact and approximated casesPl 
Then, conditioned on V, the interference Y = z^i=i ^.t 
becomes a sum of independent CS Gaussian r.v.'s and is 
therefore a CS Gaussian rv. given b}0 



y'^K{0,2AVx) 



where A is defined as 

A- 

Furthermore, we show in 
be expressed as 

E fc2 



1=1 



R 



2b 



(11) 



(12) 



that after some algebra, Vx can 



T6f = y + Y^{a^a', cos{9, - e^)}, i > 1, (13) 

where E = fc^E{af} is the average symbol energy of each 
interfering node, measured 1 m away from the interferer|3 
Because the rv.'s {X;}^]^ are i.i.d., Vx does not depend on i 
and is only a function of the interferers' signal constellation. 
For the case of equiprobable symbols and a constellation that 
is symmetric with respect to the origin of the IQ-plan^H 
(e.g., il/-PSK and 1/-QAM), the second right-hand term in 
(O vanishes and Vx ~ EjZ. 

Lastly, note that since A in ( fT2b depends on the interferer 
positions V (i.e., {Ri\°^x {Gi}^i), it can be seen as 
a rv. whose value is different for each realization of V. 
Furthermore, Appendix|B shows that r.v. A has a skewed stable 
distribution 1251 given b) 3 



A^S 



where 6 > 1, and Ct is defined as 



f~i A 



V{2-x) ca&^-KxITj ' 



x-^ 1, 



(14) 



(15) 



This distribution is plotted in Fig. [3] for different 6 and A. 



'^We can obtain (TO) following another approach: if we consider that the 
interfering nodes are coded and operating close to capacity, then the signal 
transmitted by each interferer is Gaussian, such that Xi ~ A/'c(0, 2Vx) 1231 . 

n 1^' 

We use X ~ to denote the distribution of rv. X conditional on Y. 

'''Unless otherwise stated, we will simply refer to E as the "average symbol 
energy" of the interferers. 

'^A constellation is said to be symmetric with respect to the origin if for 
every constellation point {x, y) £ "S?, the point (—x, —y) also belongs to the 
constellation. 

'^We use iS(o, 7) to denote a real stable distribution with characteristic 
exponent cx £ (0, 2], skewiiess f} G [—1, 1], and dispersion 7 6 [0, 00). The 
corresponding characteristic function is 



<j){w) ■■ 



exp [-7|«)|" (1 - j/3sign(«;)tan ^)] , a I, 
exp [— 7|tt;| (1 + j^l3sign{w) In |«>|)] , a = 1. 
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Figure 3. P.d.f. of A for different amplitude loss exponents b and interferer 
densities A (a^g = 10). Stable laws are a direct generalization of Gaussian 
distiibutions, and include other densities with heavier (algebraic) tails. 



C. Unconditional Interference Distribution 

To derive the unconditional distributiorPl of the aggregate 
interference Y in (|7]i-(l8]l, we can show that sums of the form of 
(O belong to the class of symmetric stable distributions 1251 . 
This is because the r.v.'s {Ri}^i correspond to distances in 
a spatial Poisson process and the {X;}^]^ are i.i.d. and have 
a CS distribution. Specifically, Appendix O shows that Y has 
a CS complex stable distribution given b\l"^l 



Y -5c 



7Y = XnC-y'^'/'^EilX,,^]'/"}) , (16) 

where b > 1, and Cx is defined in ( fTsT l. Using (|4li-(|5]i, we can 
further express E{|X,.„p/^} in as 



xE 



a^^cos{e, + (|)^) + a'Jl - ^ ] cos{9'i + 



2/6" 



b 


E{|X, 
E 

BPSK 


^\ ' } 
QPSK 


1.5 


0.374 


0.385 


2 


0.423 


0.441 


3 


0.509 


0.531 


4 


0.576 


0.599 



= x{b) 



(17) 

For the particular case of Rayleigh fading, ([IT) reduces to 
E{|X,,„|2/*'} ^ fc2/''r(l + i) • x{b), where we have used the 
moment relation for the Rayleigh rv.'s ||26l . Since different 
interferers i transmit asynchronously and independently, the 
parameter x(6) does not depend on i and is only a function 
of the amplitude loss exponent b and the interferers' signal 
constellation. Table J] provides some numerical values for 

n\X^,n\'/'}■ 



Table I 

E{|^i,nP''''} FOR VARIOUS AMPLITUDE LOSS EXPONENTS b AND 
MODULATIONS, ASSUMING RAYLEIGH FADING. NOTE THAT FOR M-PSK 
MODULATIONS, THIS QUANTITY IS PROPORTIONAL TO E'^/'', WHERE E IS 
THE AVERAGE SYMBOL ENERGY OF THE INTERFERERS. 



D. Discussion 

The results of this section have to be interpreted carefully, 
because of the different types of conditioning involved. In the 
unconditional case, we let V be random, i.e., we let {Ri}°Zi 
be the random outcomes of an underlying spatial Poisson 
process, and {Gi}°^i be the random shadowing affecting each 
interferer. Then, the unconditional interference Y is exactly 
stable-distributed and given by ( fT6b . We note that ( fTSI l and 
il% hold for a broad class of fading distributions, in addition 
to Rayleigh fading. In the -conditioned case, the positions 
of the interferers are fixed. Then, A in ( fT2] i is also a fixed 
number, and the interference Y is approximately CS Gaussian 
with total variance 2AVx, as given in (fTTT i. 

IV. Error Probability 

In the previous section, we determined the statistical dis- 
tribution of the aggregate interference at the output of a 
conventional linear receiver We now use such result to directly 
characterize of the error probability of the probe link, when 
subject to both interference and thermal noise, in both cases 
of slow and fast- varying V. 

A. Slow-varying Interferer Positions V 

In the quasi-static scenario of slow-varying V, it is insightful 
to analyze the error probability conditioned on a given real- 
ization V of the distances {Ri}°°^^ and shadowing {Gi}f^i 
associated with the interferers, as well as on the shadowing Go 
of the probe transmitter We denote this conditional symbol 
error probability by Pe(Go, V)^ 

To derive the conditional error probability, we employ the 
results of Section IIII-BI for the P-conditioned distribution of 
the aggregate interference Y. Specifically, using (|9]l and (fTTT i. 
the received signal Z in (|6]l can be rewritten as 



-aoe- 



W, 



where 



(18) 



(19) 



Unconditional in the sense of being averaged over the positions V. 
'^We use cSc(a, /3 = 0,7) to denote a CS complex stable distribution 
with characteristic exponent a and dispersion 7, and whose characteristic 
function is </>(w) = cxp(— 7|w|°'). Furthermore, the corresponding real and 
imaginary components are both S(a, 13 = 0, 7). 



W = Y + W ^ AAe(0, 2AVx + No), 

and A was defined in ( fT2] i. Our framework has thus reduced 
the analysis to a Gaussian problem, where the combined 

"The notation Pe(X, Y) is used as a shorthand for P{error|X, Y]. 
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noise W is Gaussian when conditioned on the location of 
the interferers. 

The corresponding error probabiHty Pe(Go, V) can be found 
by taking the well-known error probability expressions for 
coherent detection of linear modulations in the presence of 
AWGN and fast fading EtJ-JSOI, but using 2AVx + Nq 
instead of No for the total noise variance. Note that this 
substitution is valid for any linear modulation, allowing the 
traditional results to be extended to include the effect of 
network interference. For the case where the probe transmitter 
employs an arbitrary signal constellation in the IQ-plane and 
the fading is Rayleigh-distributed, the conditional symbol error 
probability Pf.{Go,T') is given by 



M 



k=l l£Bk 



1 

2^ 



WkA 



VA 



where 



VA 



„2bi 



Nn 



dO, (20) 



(21) 



'o\2AVx 

is the received signal-to-interference-plus-noise ratio (SINR), 
averaged over the fast fading; M is the constellation size; 
{Pk}iLi are the symbol probabilities; Bk, (t)k,u Wk,u and i^kd 
are the parameters that describe the geometry of the con- 
stellation (see Fig. IDi; Eq = fc^E{aQ} is the average symbol 
energy of the probe transmitter, measured 1 m away from the 
transmitter; A is defined in (fTZt and distributed according to 
(fl4] i; and Vx is given in ( fT3T l. When the probe transmitter 
employs A/-PSK and 71/-QAM modulations with equiprobable 
symbols, ( |20l i is equivalent tcF*! 



and 



M -1 



M 



(22) 



4 1- 



■Ia 



2 2{M -1) 
3 



4' 2{M -1] 



respectively, where the integral lA{x,g) is given by 
1 ^ 



^A{x,g) 



1 



9 



:VA 



de. 



(23) 



(24) 



In the general expression given in (|20] | and ( l2Tl i. the network 
interference is accounted for by the term 2AVx, where A 
depends on the spatial distribution of the interferers and prop- 
agation characteristics of the medium, while Vx depends on 
the interferer transmission characteristics. Since 2AVx simply 
adds to No, we conclude that the effect of the interference 
on the error probability is simply an increase in the noise 
level, a fact which is intuitively satisfying. Furthermore, note 
that the modulation of the interfering nodes only affects 

^"in this paper, we implicitly assume that A/-QAM employs a square signal 
constellation with M = 2" points (n even). 




Figure 4. Typical decision region associated with symbol si. In general, 
for a constellation with signal points sj. = |sfe|e^^fc and 



k = 1 . . . M, four parameters are required to compute the error probability: 
(pk.l and ipf^ i are the angles that describe the decision region corresponding 
to Sfe (as depicted); Bk is the set consisting of the indexes for the signal points 
that share a decision boundary with (in the example, Bi = {2, 3, 4}); and 



the term Vx, while the (possibly different) modulation of 
the probe transmitter affects the type of error probability 
expression, leading to forms such as ( |22] | or (l23l l. 

In our quasi-static model, the conditional error probability in 
(I20I 1 is seen to be a function of the slow-varying user positions 
and shadowing (i.e.. Go and V). Since these quantities are 
random, the error probability itself is a r.v. Then, with some 
probability. Go and V are such that the error probability of 
the probe link is above some target p*. The system is said to 
be in outage, and the error outage probability is 



„.v{Pe{Go,V)>p*}, 



(25) 



In the case of slow-varying user positions, the error outage 
probability is a more meaningful metric than the error proba- 
bility averaged over V. 



B. Fast-varying Interferer Positions V 

In the dynamic scenario of fast- varying V, it is insightful 
to average the error probability over aU possible realizations 
of interferer positions V. We denote this average symbol error 
probability by Pe(Go). Note that we choose not to average 
out the shadowing Go affecting the probe transmitter, since 
we have assumed the probe transmitter node is immobile at 
a deterministic distance j'o from the probe receiver, and thus 
Go is slow-varying. 

To derive the average error probability, we use the decom- 
position property of stable r.v.'s fl5\, which allows Y in (fTSI l 
to be decomposed as 



Y = \/SG, 
where B and G are independent rv.'s, and 



(26) 



B r-. S [as = Pb = 1, IB = cos ^ j (27) 
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G^KiO,2VG), Vg = 262-'/" (A7rC2-/,E{|X,,„p/b})\ 

(28) 

with E{|Xi_„p/''} given in ^l}. Conditioning on the r.v. B, 
we then use (|9|l and (l26b to rewrite the aggregate received 
signal Z in (|6l) as 



where 



-aoe^'" + W, 



W = /BG + W ^ AAc(0, 2BVg + No). (29) 



Again, our framework has reduced the analysis to a Gaussian 
problem, where the combined noise W is a Gaussian r.v. 
Note that this result was derived without resorting to any 
approximations - in particular, the Gaussian approximation of 
([Tol l was not needed here. We merely used the decomposition 
property of symmetric stable r.v.'s. 

The corresponding error probability Pe{Go) can be found 
by taking the error expressions for coherent detection in the 
presence of AWGN and fast fading, then using 2BVg + Nq 
instead of A^o for the total noise variance, and lastly (unlike 
in Section ITV-Ab averaging over the rv. B. For the case where 
the probe transmitter employs an arbitrary signal constellation 
in the IQ-plane and the fading is Rayleigh-distributed, the 
average symbol error probability Pe(Go) is given by 



Pe(G, 



M ^ 




where 



rf{2BVG + No)' 



d0, (30) 



(31) 



B is distributed according to (IZTT i; Vg is given in (l28T i; and the 
other parameters have the same meaning as in Section IIV-AI 
When the probe transmitter employs A/-PSK and A/-QAM 
modulations with equiprobable symbols, ( l20l i is equivalent to 



'{Go) =Ib 



and 



4 1 



'M 



(32) 



4 1 



1 



2' 2{M - 1) 

TT 3 



/MJ "V4'2(M-1) 
respectively, where the integral lB{x,g) is given by 



1 r 

^B{x,g) = - / 

Jo 




(33) 



(34) 



C. Discussion 



Using the results derived in Sections IIV-AI and IIV-BI we 
can now analyze the dependence of the error performance 
on the density A and the average symbol energy E of the 




0.01 

interferer density A (m" 



Figure 5. INR — X curves of constant P^^,, for the case of slow-varying 
interferer positions V (BPSK, SNR = 40 dB, b = 2, ro = 1 m, (TdB = 10, 
p* = 10^^). INR is tlie interference-to-noise ratio, detined as INR = E/Nq. 
Clearly, for a fixed error performance, there is a tradeoff between the density 
and energy of the interferers: if the INR (or, equivalently, E) increases, A 
must decrease, and vice-versa, to maintain the same outage probability. 



interfering nodes. For that purpose, we use ( |20] |. although 
(l30l l would lead to similar conclusions. In ( l20l i. the error 
probability Ps{Go, V) implicitly depends on parameters A and 
E through the product AVx in the denominator of t^a in ( |2TI ). 
This is because the dispersion parameter 7^ of the stable rv. A 
depends on A according to (fT4l) . and Vx is proportional to E as 
in ( fTsT l. The dependence on A can be made evident by using the 
scaling property of stable r.v.'s [25] to write AVx — ^^AVx, 
where A is a normalized version of A, independent of A. We 
thus conclude that the interference term AVx is proportional 
to X'^E, where 6 > 1. Clearly, the error performance degrades 
faster with an increase in the density of interferers than with 
an increase in their transmitted power. The tradeoff between 
E and A for a fixed error performance is illustrated in Fig. |5] 



D. Numerical Results 

Figs. |6] and [T] quantify the average and outage proba- 
bilities for several scenarios, showing their dependence on 
various parameters involved, such as the signal-to-noise ra- 
tio SNR = Eq/No, interference-to-noise ratio INR = E/No, 
amplitude loss exponent b, interferer spatial density A, and 
link length tq. 

The plots of Pgy, and Pe{Go) presented here are of semi- 
analytical nature. Specifically, we resort to a hybrid method 
where we employ the analytical results given in (I20l)-(l25b 
and (I30]l-(l34li. and perform a Monte Carlo simulation with 
respect to the stable r.v.'s (i.e., A and B), according to ||3TI . 
Nevertheless, we emphasize that the expressions derived in 
this paper completely eliminate the need for simulation of the 
interferers' position and waveforms in the network, in order 
to obtain the error performance. 

For illustration purposes, the plots assume that all terminals 
(i.e., the probe transmitter and interfering nodes) use BPSK 
modulation. We analyze both cases of slow and fast-varying 
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10 15 20 25 30 35 40 

SNR (dB) 



(a) P^u, versus the SNR of the probe link, for various interference-to-noise 
ratios INR (BPSK, b = 2, X = 0.01 m'^, ro = 1 m, (t<jb = 10, p* = 10"^) 




10 15 20 25 30 35 40 

SNR (dB) 



(b) P^j,, versus the SNR of the probe link, for various interferer spatial 
densities A (BPSK, INR = 10 dB, 6 = 2, ro = Im, o-jb = 10, p* = 10~^). 

Figure 6. Error outage probability plots for a heterogeneous network 
(where SNR INR in general) and slow-varying interferer positions P. 
Since P is slow-varying, the meaningful performance metric is the outage 
probability P^^ given in )25t . 



interferer positions V, concurrently with the following two 
different scenarios; 

1) Heterogeneous network: The probe transmitter is al- 
lowed to use an arbitrary power Pq = Eq/T, not 
necessarily equal to the common power of the interfering 
nodes P = E/T, and hence SNR ^ INR in general. This 
scenario is useful when the goal is to evaluate the impact 
of aggregate interference from a large number of identi- 
cal secondary users (e.g., cognitive-radio terminals) on 
the performance of a primary link. 

2) Homogeneous network: The probe transmitter and in- 
terfering nodes all use the same power, and thus 
SNR = INR. This may correspond to a sensor network 
scenario, where there is a large number of indistinguish- 
able, spatially scattered nodes with similar transmission 
characteristics. In such a case, the goal is to evaluate the 
impact of the aggregate network self-interference on the 




probe link length ni (m) 



(a) Pe(Go) versus the length r'o of the probe link, for various signal 
loss exponents b (BPSK, Go = 0, SNR = INR = 20dB, A = 0.01m-2^ 
TdB = 10). 




5 10 15 20 25 30 

SNR = INR (dB) 

(b) Pc(Go) versus the SNR, for various interferer densities A (BPSK, Go = 0, 
fe = 3, ro = 1 m, (TdB = 10). 



Figure 7. Average error probability plots for a homogeneous network (where 
SNR = INR) and fast- varying interferer positions P. Since P is fast- varying, 
the meaningful performance metric is the average eiTor probability Pc(Go) 
given in i30\ . For simplicity, we use Gq = in these plots (no shadowing 
on the probe link). 

performance of each sensor node. 
For the heterogeneous case depicted in Fig. |6] we conclude 
that the error performance deteriorates as A or INR increase, 
for a fixed SNR. This is expected because as the density or 
transmitted energy of the interferers increase, the aggregate 
interference at the probe receiver becomes stronger. Note, 
however, that in the homogeneous case where SNR = INR, 
the error performance improves as we increase the common 
transmitted power P of the nodes (or equivalently, the SNR), 
although the gains become marginally small as P oo (see 
Fig. |7(b)[ ). This happens because in the interference-limited 
regime where SNR = INR ^ 1, the noise term A^o in <EB or 
(ISTT i becomes irrelevant, and so the SNR in the numerator can- 
cels with the INR in the denominator, making the performance 
independent of the transmitted power P. 

The effect of the amplitude loss exponent b on the error 
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performance, on the other hand, is non-trivial. As illustrated 
in Fig. |7(a)[ an increase in b may degrade or improve the 
performance, depending on the value of the link length j q and 
other parameters. This is because b simultaneously affects both 
the received signal of interest and the aggregate interference - 
the former, through the term I/tq; and the latter, through a a 
and in (fT4l i. or through as, 7b, and Vq in jZTI i and ( |28] l. 



V. Summary 

This paper introduces a mathematical model for commu- 
nication subject to network interference and noise. The in- 
terferers are scattered according to a spatial Poisson process, 
and are operating asynchronously in a wireless environment 
subject to path loss, shadowing, and multipath fading. We 
show that the aggregate network interference at the output 
of a linear receiver is related to a skewed stable distribution 
when conditioned on the positions of interferers, and to a 
symmetric stable distribution in the unconditional case. We 
characterize the error performance for the cases of slow and 
fast-varying interferers, in terms of outage and average error 
probabilities, respectively. These expressions are valid for any 
linear modulation scheme. We then quantify these metrics as 
a function of various important system parameters, such as the 
SNR, INR, path loss exponent, and spatial density of the inter- 
ferers. In Part II of the paper |[T3l . we characterize the capacity 
of the link when subject to both network interference and 
noise, and derive the spectrum of the aggregate interference 
at any location in the plane. Lastly, we put forth the concept 
of spectral outage probability, a new characterization of the 
aggregate interference generated by communicating nodes in 
a wireless network. 



Appendix A 

Derivation of the Complex Baseband Interference 
Representation 

To derive the representation ^ and (O of the aggregate 
interference Y{t), we project Y{t) onto the basis func- 
tion Tpiit) = ^2/rcos(27r/ct) as follows: 

Yi= [ Y{t)i>i{t)dt 



R 



—at cos(27r/ct + 6t + - t) 



+ \l -a', cos(27r/et + 91 + (t,,)u{t - A) 



X -\/ — cos(27r/ci)dt 



^ 2 ka.e"^^ 



2 Jo 



Di 



cos{9i + (j)i)dt 



Oh 
2" 



Di 



I rT 



cos(47r/ci + 9i + (j)i)dt - 



+ ^ / cos{A^ Ut + 9[ + (j},)dt 

^ J Di 



E 



«0 for /cT»l 
=<yGi V 



R° 



where 



a,^ cos{9, + + a\(^l- cos(0^ 



The signal Y{t) can be projected onto the basis func- 
tion ^2it) = — ■\/2/rsin(27r/c<) in an entirely analogous 
way, leading to 



where 

2 = kai 



Y. 



E 



1 - — I sin( 



We can combine Xii and Xi2 in the complex 

rv. X, = Xi^i +jXi^2 as 



ka,e3'^' 



T 



a.eJ"- + 1 



A 

T 



which completes the derivation. 

Appendix B 
Derivation of the Distribution of A 

To derive the distribution of A given in ( fl4] i. we start with 
the following theorem. 

Theorem B.l: Let {Ti}°^^ denote the arrival times of a 
one-dimensional Poisson process with rate A; let {VFjj^i 
be a sequence of nonnegative i.i.d. rv.'s, independent of the 
sequence {r^} and satisfying E{|VFi|"} < oo. If < a < 1, 
then 

£ 5 («, /3 1, 7 = xc-^n\w.n) , 

1=1 '''i 

where Cq. is defined in ( fTSl ). 

Proof: See EH. □ 
If an homogeneous Poisson point process in the plane has 
spatial density A, and Ri denotes the distance of node i 
to the origin, then the sequence {-Rfj^i represents Poisson 
arrival times on the line with the constant arrival rate Att. This 
can be easily shown by mapping the spatial Poisson process 
from Cartesian into polar coordinates, and then applying the 
mapping theorem 1T4| . Using this fact, we can then apply the 
above theorem to ( fT2b and write 



X ^ e 



2(jGi 



i—l ^ 




cos{9- + (f>i)dt 



Di 



5 fa = i /3 = 1, 7 = ATrC-iElle^-^-li/"}) , (35) 



bOfor/,T»l 
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for b > I. Using the moment property of log-normal rv.'s, i.e., 
E{e''^} = e'='/2 for G - 7V(0, 1), ^ simplifies to 

for 6 > 1. This is the result in (fT4l i and the derivation is 
complete. 

Appendix C 
Derivation of the Distribution of Y 

To derive the distribution of Y given in (fT6] l. we start with 
the following theorem. 

Theorem C.l: Let {Ti}°^i denote the arrival times of a one- 
dimensional Poisson process with rate A; let {Zi}^-^ be a 
sequence of CS i.i.d. complex r.v.'s Z,; = Zi^i + jZi^2, inde- 
pendent of the sequence {ri} and satisfying E{|Zi|"} < oo. 
If < a < 2, then 

E ^ (a, /? = 0, 7 = XC-'E{\Z,,r}) , 

1=1 "^i 

where Cq, is defined in ( fT5] ). 

Proof: See 1251 . For an alternative proof based on the 
influence function method, see 1321 . □ 
Using the Poisson mapping theorem as in Appendix |BJ we 
can apply the above theorem to (|7]i and write 



E 

i=l 



E 




0, 7 = A7rC27,E{|e^^'X,,„|2/!>}^ 



(36) 

for b > 1. Note that Xi, whose expression is given in (O, 
is CS due to the uniform phase (pi. As a result, e'^'^^'X^ is 
also CS. Using the moment property of log-normal rv.'s, i.e., 
E{e''^} = e''^'/^ with G - A^CO, 1), dMI simplifies to 



Y - 5c 



(« = ^, /3 = 0, 7 = A^C-/,e2<xVb^E{|X,,„|2/''}^ , 

for 6 > 1. This is the result in (fTST i and the derivation is 
complete. 
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